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We investigate the stationary nonequilibrium states of a quasi one-dimensional system 
of heavy particles whose interaction is mediated by purely elastic collisions with light 
particles, in contact at the boundary with two heat baths with fixed temperatures 
, T + and T~ . It is shown that Fourier law is satisfied with a thermal conductivity 

■ proportional to y/T{x) where T(x) is the local temperature. Entropy flux and entropy 

production are also investigated. 
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I. INTRODUCTION 



Although the Fourier's law ]q = — k gradT, relating the macroscopic heat flux jq to the temper- 
ature gradient gradT, has been introduced almost two centuries ago [1], its microscopic basis is still 
an open issue. Actually, its status has long been purely phenomenological and its only justification 
was then its accurate and faithful agreement with observations in numerous experimental instances. 
Following the advances of nonequilibrium statistical mechanics, several microscopic models has been 
recently introduced aiming at deducing the Fourier's law from microscopic principles and relating the 
conductivity to the microscopic parameters of the system. In particular, the main issue is to determine 
what are the minimal requirements to get Fourier's law [2] [3] [4] [5] [6] [7]. 

The present study belongs to this line of research. It exploits recent results [9] [10] about the 
so-called "adiabatic piston problem" to device a spatially extended, quasi-one-dimensional system, 
in contact at its ends with two heat baths at different temperatures (T~ on the left and T + on the 
right) . It is composed of an array of K "pistons" (heavy point particles) separated by compartments 
filled with A7V non-interacting light particles of mass m. The actual interactions are limited to 
elastic collisions between light particles and adjacent pistons: successive pistons are thus coupled 
through their interactions with the same fluid, while at the same time successive fluid compartments 
are coupled through their interactions with the same intermediary piston. The differing boundary 
temperatures force the system out of equilibrium, and one of the issues tackled in this paper is to 
determine the temperature profile in the system. A key point, discussed in [8] and [10], is that as 
soon as temperatures differ on each side of a piston, its stochastic motion induces an heat transport 
(whereas the pistons are adiabatic when fixed) . Our model is simple enough to determine a consistent 
stationary state (with no drift) for any given heat flux; it thus sheds light on the controversial issue 
of heat conductivity in 1-dimension [11]. The derivation is performed analytically, in the frame of a 
perturbation approach developed in [9] with small parameter m/AM where m is the mass of fluid 
particles and AM the mass of one piston. Such an array is somehow reminiscent of chains of masses 
linked by springs, much studied since the pioneering work of Fermi, Pasta and Ulam [12] [13]. The 
masses correspond here to the pistons and the springs to the fluid compartments. But rather than 
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being the nonlincarities, the key ingredient at the origin of transport in our setting will appear to be 
noncquilibrium fluctuations. Our model can also be related to the general class of models introduced 
by Eckmann and Young [7] with the pistons playing the role of the "Energy Storing Devices", and 
the fluid particles the role of the "Conducting Agents" . 

The paper is organized as follows. In Section II, we describe the stationary state and heat flux for 
a single piston surrounded with two fluids at different temperatures. In Section III, we investigate 
the fluid separating two pistons and the indirect coupling of the pistons that it achieves. In Section 
IV, we then bridge these two sets of results to get the behavior and heat conductivity of the above- 
mentioned array. In Section V, the Fourier's law is recovered and discussed. In Section VI, we give 
explicit relations for entropy production and dissipation, before ending the paper with a final Section 
VII devoted to a discussion of the scope of our results. 

II. SINGLE PISTON STATIONARY STATE WITH HEAT FLUX (AND NO DRIFT) 

In [10] we investigated the quasi 1-dimcnsional problem of non-interacting point particles of mass 
m colliding elastically with a single heavy "piston" of mass M ^> m. Initially the piston is fixed at 
X = and the particles on the left (resp. right) of the piston are in thermal equilibrium described 
by a Maxwcllian distribution of velocities (parallel to the axis of the cylinder*) with temperature T ~ 
and uniform density tIq , i.e. pressure p^ — UqUbTq (resp. Tq~, Uq , Pq). At time trj the piston is let 
free to move without any friction. The initial conditions (pg , T^) were chosen such that the piston, 
which moves stochastically under the collisions with the particles, remains on the average at X = 0. 
It was then shown that the system evolves to a noncquilibrium stationary state (Fig. 1) where the 
piston has a temperature Tp (average kinetic energy) and the fluids on the left/right of the piston are 
characterized by temperatures T~ ^ T+, pressure p~ = p + = p, heat current ]q = jg = jq, and no 
drift (w~ = w + = 0). 
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Figure 1: Stationary state of the piston problem (with no drift). If Tq / T , it is a nonequilibrium state with 

jq + and T+^T". 



A. Stationary state of the fluid with prescribed {p,T,jq) and no drift 
It was furthermore shown in [10] that the nonequilibrium stationary state of the fluid parametrized 

by 

(p, T, j Q ,w = 0) (1) 
is characterized by the bimodal distribution function (9 being the Heaviside function) 



"The other velocity components, playing no role in this problem, can be initially set to zero and then remain 
zero at all times. 
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p{x,v,t) = P (v) = 2 J> JH-fJ* e-^ + 9[ v] ——- \ j~ e-^/v (2) 

m I 1 + y V 7r 1 + y y 7ry J 

where 

m , % 

^2^T < 3 > 

(not to be confused with the inverse temperature l/ksT: here /? -1 / 2 is a velocity) and 

y = i + ^-fV4 + c^ (4) 

with C the dimensionless parameter 

C = V^/3 1/2 ^ (5) 
P 

Let us note that the parameter y is strictly positive. Moreover since 

then j'q = —jq implies y' = 1/y. It is straightforward to check the consistency between the above 
result (2) for the velocity distribution and the physical meaning of the associated thermodynamic 
parameters: 

/dv p(v) = n = (number density n of the ideal fluid) (7) 

k B T 

J dv p(v)v = (no drift) (8) 

dv p(v)v 2 — (kinetic temperature T of the fluid) (9) 



/ 



2 

dv p(v)v — — jq (heat flux j'q) (10) 



m 



It is to be underlined that Eq. (2) reflects the fact that three different temperatures are involved in the 
nonequilibrium stationary state of the fluid (see Fig. 1): the temperature T ~ = T/y of particles going 
to the right (v > 0), the temperature = Ty of particles going to the left (v < 0) and the average 
kinetic temperature T defined by (9). These three temperatures coincide only at thermal equilibrium, 
when y = I, jq =0 and p(v) coincide with the usual Maxwellian distribution. They differ as soon 
as a heat current jq ^ is forced into the system and shifts the velocity distribution away from the 
equilibrium Maxwellian distribution. Note that y < 1 iff jq is positive, and y > 1 iff jq is negative. 
The parameters of the fluid in the nonequilibrium state can thus be reformulated as 

(p, p, y = y{p,(3,j Q ), w = 0) (11) 



B. Stationary state of the piston 

In the situation sketched in Fig. 1, we denote by T~ (rcsp. T+) the temperature of the fluid in 
the stationary state on the left (resp. right) side of the piston, and by Tp the temperature of the 
piston. In a stationary state, Eq. (2) is to be written on each side of the piston with respectively 
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(3 = j3~ = m/2kBT~ and y~ = y(p, f3~ , ]q) in the left compartment, and (3 — (3 + — m/2kBT + and 
y + = y{p 1 jq) in the right compartment. Xp being the position of the piston, it was shown in [10] 
that in the stationary state for the piston 



T(X P - 0) = TqVI + 5- (left of the piston) 
T(X P + 0) = T+Vl + 5+ (right of the piston) 



(12) 
(13) 



where 



(14) 



with Tp being the temperature of the piston (average kinetic energy of its stochastic motion) and 

2m 



a = 



M + m 



(15) 



C. Stationarity state of the piston and surrounding fluids 

Stationarity for the whole system (piston, fluid on the left and fluid on the right) finally yields the 
consistency conditions 



T(X P - 0) = T~ = T ~y 
From (4) we have 



T+ 

and T(X P + 0) = T+ = -2- 

V 



■--=-C V4 + C 2 



(16) 



(17) 



and thus from (13) and (14) 



T P = T+ 



y+ + 



a(2 — a) \y+ 



-y 



Similarly, from (14) 

T P = T 

= T 



1 1 

y- a(2 - a) \ V y~ 



- 1 



1 + ^» 2 - C V 1 + ^Va(2- ( ,, 



- 1 



(18) 



(19) 



From the definition of C, Eq. (5), and the fact that p + = p = p and jg = jg = jq, we have the 
equality T + (C + ) 2 = T~{C~) 2 relating the fluid nonequilibrium parameters on each side of the piston 
in the stationary state; this relation gives, using (18) and (19) 



/ T + 



l T - = _ 



a{2 - a) 



(20) 
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with 



2 , _ M 

a(2 — a) 2m 



1 



M 2 



(21) 



which shows that jq = 0(m/M), as discussed in [10]. Let us remark that (18) and (19) imply also 
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2T P = T + +T~ + 



JQ 
2k B V P 



+ 



+ 



a(2 -a) 2 



7TTO 

2k B p 



f T+\U + 



Trm 



JQ 

2k B T+ V P 



-VT-4 4 + 



2k B T- 



(22) 



Using Eq. (20) and the fact that jq = 0(m/M), we conclude that the temperature of the piston 
separating the two fluids is given by 



T P = VT+T- + O 



"/my 



(23) 



III. STATIONARY STATE OF TWO OR MORE PISTONS 

Let us now consider the stationary state defined by (p, T, Jq,w — 0) where the fluid is bounded 
by two identical stochastic pistons, each of mass M 3> m, which remain on the average at the same 
position under collisions from both sides (no drift). We shall use the same notations as in [10] where 
Tq (respectively Tq~) denotes the temperature of the particles incident on the right piston from the 
left, i.e. with v > (respectively on the left piston from the right, i.e. with v < 0), and T is the 
(average kinetic) temperature of the fluid, Eq. (9), in the intermediary compartment in the stationary 
state, see Fig. 2. 
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Figure 2: Stationary state with two pistons with prescribed current jq, prescribed boundary temperatures T , 
and no drift (uniform pressure p). T ~ and Tq are auxiliary temperatures, associated with subpopulations of fluid 
particles, going respectively to the right (v > 0) and to the left (v < 0) inside the middle fluid compartment. 



A. Temperature Tp of right piston 



From Eqs. (12) and (16), we have on the left side of the right piston 

T T 

y 



T ~ = - 



VT+6 

(joint stationarity of the fluid and the right piston) i.e. 



(24) 
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V 2 = 1 + 6- 



(25) 



Denoting Tp the temperature of the right piston (see Fig. 2), plugging Eqs. (24-25) in Eq. (14) for 
5~ yields 



y 2 - 1 = a(2 - a) 



(26) 



i.e. 



1 



y a(2 - a) 



(27) 



B. Temperature T P of left piston 

Similarly denoting Tp the temperature of the left piston (see Fig. 2), Eqs. (13) and (16) give for 
the right side of the left piston 



T + = Ty 



VTT5+ 



(28) 



i.e. 



+ ^2 



- 1 



(29) 



and thus from Eqs. (14) and (28-29) we obtain 



1 



a(2 - a) \y 



- -y 



Let us note that 



7 /' 4 T/: [ [ - ■ n = 77:2 + C 1 ) 



(30) 



(31) 



Therefore the temperature T of the fluid is related to the temperature of the surrounding pistons by 
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T-^Tp+Tp) —y- 



(32) 



Results of Section II can then be applied to each piston to obtain the temperature T~ on the left of 
the left piston and temperature T + on the right of the right piston. Bridging the ensuing relations 
allows to determine heat conductivity of this composite system, as follows in the next section. 



IV. EXPRESSION OF HEAT CONDUCTIVITY 



A. Elementary case of a single fluid compartment bounded by two pistons 



We first derive auxiliary relations for an elementary unit, by considering again the stationary state 
for two pistons of area A separated (on the average) by a distance AL, with A/TV points particles 
between them (see Fig. 2). From (27) and (30), we have 
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Tj?-T£ = T 



1 



y 



1 



a(2 — a] 



(33) 



Introducing Eq.(17) and the expression (5) of C, together with (21) and p A AL = AN ks T, we 
obtain 



M ( m 2 \ nun ( A \ 1 / 
m[ 1+ M^J\l2k% jQ [ANj jfV 



\AN) VT 



8k B T 



where the temperature T of the fluid between the two pistons is given by (32), i.e. 

T--(T P +T P ) —y- 
For the fluids on the left and on the right, Eq. (23) gives at lower order in m/M 
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T+ = + O 



/my 



(34) 
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Figure 3: Temperature profile for the stationary state of two identical pistons. Tp and Tp are the temperatures 
(average kinetic energy) of the pistons, whereas T~, T and T + are temperatures within the fluids as defined in (9). 



B. Heat transport through an array of K pistons 

We are now in position to derive the law describing the fluctuation-driven heat transport equation 
for our model of spatially extended system and to investigate whether, and in which conditions, the 
usual Fourier's law can be recovered. We now consider the system of K identical (movable) pistons 
of mass AM, with AN point particles of mass m between them. In the stationary state defined by 
(T~ ,p,jo) on the left, (T + ,P,Jq) on the right, the pressure and the heat flux in each compartment 
(defined by successive pistons) will also be given by p and ]q, and we denote respectively 7\ and Tp^ 
the temperatures of the successive compartments and successive pistons (see Fig. 4). The average 
distance AL k between adjacent pistons will adapt to ensure that the pressure is actually homogeneous 
across the system (local mechanical equilibrium); it is thus given by 



{Xk+i) — {Xk) 



p,i Tpk Tp,k+i TpK 
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Figure 4: Stationary state with K pistons. The pressure p and the heat flux jq are identical in each subsystem, 
the drift velocity is zero and the pistons remain on the average at fixed positions (Xk). 



From Eqs. (34-35) applied to the fc-th compartment 



1 



AL k 
with 



P,k+1 



Tp,k) 



AM 



m 



1 + 



' m y 



7rm / A \ 
~AN J 



2kg 



JQ 



'1 + 



nm 



Sk B T k \ p 



^ ) (38) 



T k = o( T P,fe + Tp.k+i) - ^— ( 



(39) 



The system considered is defined by the total mass of pistons M = KAM, the total number of point 
particles N = (K — 1)AN, the uniform pressure p and the temperatures T~ and T+ of the outer 
(semi-infinite) compartments. With K sufficiently large, but such that 



AM = — > m, 



1 < AN 



N 



K-l 



ALp 



(40) 



Eqs. (38-39) give 



dTp(x) 
dx 



M A 
rn N 



™ JQ_ 

2k% VT 



8k B T k 



(41) 



T(x)=T P {x) 



1 - 



4,k B T F 



where the position of the piston and the location of nearby fluid can be identified, and the discrete 
set of temperature values interpolated by smooth functions Tp(x) and T(x) at this spatial scale and 
perturbation order. With N/A given, then jq will be of the order m/M in order that grad T is 
uniformly bounded. Therefore to first order in m/M, we have 



' dTp(x) 
dx 



k{x] 



JQ 



(42) 



k(x) — 



m N 



M A V irm 



2k% \/W 



which gives by integration the temperature profile in the system 

Tp(x) - T(x) = [(T-f 2 + | ([T+f 2 - (T-) 3 / 2 ) 



2/3 



This temperature profile fulfills, in a purely Hamiltonian setting, the expected general relation 

t(x) = \(T-) a + y ((T+r - (T-ry 



(43) 



(44) 



given in [7], hence supports its universality. We here find an exponent a — 2/3 as mcntionncd in [7], 
related to the fact that the fluids are confined to a limited region between pistons and the energy 
is purely kinetic. Our result moreover shows that local thermodynamic equilibrium is not required 
to get such a temperature profile provided there is a natural way to define local temperature, here 
through (2), see also [10]. We recover the currently observed linear profile if \T + — T~\ <C T~, or in 



8 



boundary layers where x <C L or (L — x) <C L (there interchanging the role of T + and T in (43)). 
From Eq. (42), we finally get 



dT m N 2k% 2 



JQ K dx M A \l nm 3L 



(T+f/ 2 - {T-f 2 } (45) 



Moreover from the knowledge of the integrated number density profile N(x) of the fluid, such that 
JV(0) = 0, N(L) = N and (37) 



dN Ap Ap 



we obtain the length L of the system: 



Nk B = SpAL {T+ ^ 2 J_~ )3/2 (47) 



p A L = [T+ + T~ + Vf+T^} (48) 



Therefore 



^) = -^V — PWT*-^] (49) 
M V nm 

which shows that the term we have neglected when passing from (41) to (42) is of the order (m/M) 2 . 

V. FOURIER'S LAW 

We have thus established that for our model, the Fourier's law 

j Q = -K grad T (50) 

is valid at some mesoscopic scale allowing to consider the continuous limit of the discrete microscopic 
temperature profile, with the heat conductivity n given by 



M A V 7rm y ' 

We thus find that k depends on the local temperature in agreement with the general statement 
made in [7] when communicating agents (here the light particles) remain confined (here between the 
adjacent pistons). As one could expect, this heat conductivity is proportional to the number N/A 
of gas particles per unit area and also proportional to the ratio m/M, measuring the efficiency of 
interaction between the piston and gas particles. It is however independent of the length L of the 
system (except for L-dependent corrections to the continuous limit involved in the derivation of the 
gradient term). Let us note that for a strictly 1-dimensional system (for which the area A is zero), 
one should introduce the power of heat Pq which is transmitted through the system. In this case, 
Fourier's law reads 

Pq = — k grad T (52) 

and we have 



Finally, let us recall that Clausius-Maxwell-Boltzman obtained a theoretical expression for n for gases 
with k <~ VT independent of the gas density [6] . k ~ y/T for the conductivity [6] . 
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VI. ENTROPY PRODUCTION AND DISSIPATION 



Since the system considered (defined by the (K — 1) elementary units and K pistons) is in a 
stationary state, its thermodynamic entropy S remains constant. Thus, in the framework of irreversible 
thermodynamics, we have 



dS 



— -I + PqI — - — 



1 



1 



T+ 







which gives for the total entropy production per unit time of the system 

Nm 



T+ 



M 



2fc| JL ( T +_ r -)2 T+ + T-+VT+T- 
nm 3L 



T+T- 



/T+-VT- 



(54) 



(55) 



Therefore I > as required by the Second Law of thermodynamics. The same argument remains 
valid locally where the entropy density s(x,t) satisfies 



d 



d 



dt 



dx 



^-s(x,t)=i(x,t)-i-( 3 S 



(56) 



(since there is no drift, hence no extra entropy current coming from transport of matter). For the 
stationary state, it implies [10] for the source term i(x) associated with irreversibility 



i(x) 



d_ 

dx 



3Q_ 
T 



(57) 



i.e. (recall that jq is constant throughout the system) 



i{x) 



- 3Q 7^2 § rad T 

«(^) 2 



Nm 1 



M A\ m 



A(( T+ )3/2_ (T - ) 3/2 ) 



y/T(x) 



(58) 



Let us note that I is also the total entropy production for the system plus the two "reservoirs". In 
conclusion, the irreversibility i(x) is proportional to the parameter R = Nm/M. Moreover, we see 
that if R is large, which corresponds to the strong damping case [9] [14], i.e. a high dissipation, the 
heat conductivity is large whereas for R small, which corresponds to the weak damping regime, i.e. a 
weak dissipation, the heat conductivity will be negligible. 



VII. CONCLUSIONS 



We have considered a quasi-one-dimensional system composed of K — 1 (f>l) identical subsys- 
tems, in contact at the boundaries with two fluids with temperatures T^ 1 and the same pressure p. 
We have shown that there exists a nonequilibrium stationary state with heat current jq ^ given by 
(49) for which the Fourier's law jq = — k gradT is satisfied. For this model, the heat conductivity k 
is proportional to VT, as could be expected from kinetic theory. It is also proportional to the factor 
R = Nm/M, with Nm the total mass of the light particles and M the total mass of the heavy parti- 
cles (or pistons), as expected from previous work where it was shown that dissipation (i.e. friction) is 
proportional to R. It should be remarked that n is of the order of the small parameter e = m/M be- 
cause the length unit introduced is microscopic. In this stationary state, the fluid in each subsystem is 
characterized by the thermodynamic parameters (T(x),p,jQ,w — 0) where T(x) is defined as the local 
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average kinetic energy, see Eq. (9). However, it is not a state of local equilibrium since the velocity 
distribution function is given by the bimodal expression Eq. (2) and not by the Maxwell distribution 
with temperature T{x). Relaxation from an arbitrary initial condition towards this stationary state 
with no drift has not yet been investigated, but is expected to hold as was the case for one piston. The 
fluids arc moreover characterized by the equation of state p = nfe^T, e = p/2. Furthermore, the total 
length L of the system will adjust to be given by p A L = (Nk B /3) {T+ + T~ + y/T +T~). Therefore 
introducing an "average temperature" for the whole system T Av = (1/3) (T+ +T~ + VT+T"), we have 
the equation of state for a single component general ideal gas [15] pV — NUbTav, E = NksTAv/^ 
(since the number of fluid particles is much large than the number of pistons). Finally we have shown 
that the total entropy production per unit time I of our system (i.e. K — 1 subsystems) as well as the 
total irreversibility i(x) are strictly positive in agreement with the Second Law of thermodynamics. 
As mentioned in [10] , the entropy of the "heat baths" will change only due to the heat flow and there 
is no internal entropy production. Therefore the entropy production for the total system (i.e. system 
+ heat baths) is simply the quantity / computed above. 

We have assumed that the pistons are much heavier than the fluid particles. One may wonder 
whether our results are more general and remain valid if this assumption is dropped. It is an open 
problem, but clearly the mass of the "piston" must be different from the mass of the fluid particles, 
otherwise we would loose the stochasticity necessary for dissipation. It is to note that the respective 
roles of fluid compartments and pistons are more symmetric than it might seem: two successive fluid 
compartments are coupled through the fluctuations of the intermediate piston, and two successive 
pistons arc coupled through their interaction (during collisions) with the same intermediate fluid 
compartment (also fluctuating and out of equilibrium). We have here considered ideal fluids, i.e. non- 
interacting light particles. In the case of a single piston, simulations with non-interacting particles 
or with hard-core particles gave (surprinsingly) similar results [14] [9], allowing us to conjecture that 
the results given in this paper would remain valid in the case of hard-core fluid particles. Behavior 
for particles interacting through a binary finite-range potential, for which energy no longer reduces to 
kinetic energy, is still unknown [7]. 

We should stress that we have not limited our discussion to states which are close to equilibrium 
(T + — T~ is arbitrary). No assumption of linear response has thus been involved in the derivation 
of Fourier's law. In agreement with a caveat first underlined by Van Kampen [16], macroscopic 
linearity of the response is not the reflection of a linearity of the microscopic response, but rather an 
emerging feature following from averaging and cancellation of nonlinear microscopic contributions. In 
the present case, it is the perturbation approach and the underlying scale separation m/M < 1 which 
leads to cancellation of nonlinear terms at lowest order. 
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